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Abstract 

In this work we calculate the lifetime of quarkonium moving with velocity 
V through a quark gluon plasma at temperature T. We also investigate the 
I stability of heavy mesons with respect to the effects of color charge screening. 

' An explicit, configuration-space potential is found for the screened interaction 

■ between the quarks constituting the meson. We solve the Schrodinger equa- 
I tion for the relative motion of the quarks in this non-spherical potential. In 

this way, we determine the range of v, T values for which the meson is bound. 
When a bound state exists, we use the bound-state wavefunction as the initial 
state for the dissociation of the meson due to gluon absorption. The meson 
lifetime is thus determined as a function of v and T, and conclusions are 
' drawn concerning the possibility of detection of the meson in a high-energy 

heavy-nucleus collision. 

■ I. Introduction 

Calculating the lifetime of a heavy meson in a medium of quarks and gluons 
is important for the understanding of the quark-gluon plasma. If a heavy-nucleus 
collision were to produce fewer than the number of mesons expected, we would need 
to ask about the lifetime of the meson in comparison with the lifetime of the quark- 
gluon plasma. If the lifetime of the quark-gluon plasma is greater than the lifetime 
of the meson in the quark-gluon plasma, one would think that there was probably a 
suppression of mesons due to the effect of screening or due to the effect of collision 
with gluons. In fact, many claim that suppression of J/ip in heavy-ion collision 
could be a signature of quark-gluon plasma. 

There is a relationship between the inverse screening length, niei] the number 
of flavors in the quark-gluon plasma, Nf] and the temperature of the quark-gluon 
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plasma, T. If we increase the temperature of the quark-ghion plasma, the quarks and 
the gluons will become more active, and more effective at screening the interaction 
between a quark-antiquark pair. This relationship can be algebraically expressed [1] 
as follows: 

ml-\9\N + ^)T\ (1) 

where A'" = 3 from the color SU(A^ = 3) group, g is the dimensional coupling 
constant of the field strength, F^'' = d^'A"^ - d'^A^ - gfabcA'^A''^, and A^/ = 3 is the 
number of light flavors in the quark-gluon plasma, which are the up, the down, and 
the strange quarks. We will also use the temperature-dependent running coupling 
constant of QCD, which is given by [1, 2] 
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This equation exphcitly displays asymptotic freedom: — > as T — > oo. We notice 
that there is no intrinsic coupling "constant" on the right hand side of this equation. 
The only free parameter of the theory is the QCD energy scale. A, whose numerical 
value is dependent on the gauge and on the renormalization scheme chosen. If we 
choose the QCD energy scale to be A = 50 MeV [1], Nf = 3, and A^ = 3 we get 

mei = / = • (3) 

^31og(T/50) 

II. Screened Potential 

Chu and Matsui [3] studied the effects of Debye screening of the interaction 
between a heavy quark-antiquark pair moving through a quark-gluon plasma. For 
our purposes, their results may be summarized in the form of an effective in k space 
quark-antiquark screened interaction: 

where the quantities (, A, B,C, D are defined in terms of cos = x by 
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V is measured in units of c, and k is measured in units of nieic/h, the inverse Comp- 
ton wavelength associated with the screening mass, niei- The resulting interaction 
potential is given in units of ^agmeic^, with determined by QCD sum rules and 
lattice QCD calculations to be 0.232 [5]. The quantities (, A, B,C, D defined here 
are non-negative in the integration region. V(k) is invariant under rotation about 
z, and under reflection across the x — y plane. 

We obtain the quark-antiquark r-space interaction potential, U{r), by Fourier 
transforming V(k): 

f/(r) = ^Jd''k e*^-"- Vik) (5) 

In order to solve the quark-antiquark Schrodinger equation in r space, we need an 
expansion of the interaction potential in the form 

U{r)=J2ue{r)Pe{coser) (6) 

e 

If we substitute the multipole expansion of the plane wave 
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= ^7^Y.^'3,{kr)Y.{Y',{k))Xi^) 

into Equation (jSJ, the axial symmetry of ^(k) implies that only the = terms 
will survive the 0k integration. Comparison with Equation © yields 

M^(r) = /-(2£ + 1) / dxPiix) k'^dkje{kr)V{k,x). (7) 
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The integration variable x in Equation ((7j) represents cos6'k. Because V{k,x), given 
by Equation Q, is an even function of x, only even values of £ will occur in the 
multipole expansion. 

The fc-integration in Equation ((Tj) can be done exactly, using the theory of 
residues. We start with explicit expressions for the even-£ spherical Bessel functions: 

, , ^ (i + m-l)\ ai^ sm(kr) + bl, cos(kr) 

ji{kr) = 2^ — ^ — TTTz 7T — X 

m=l,2,...,e+l 

with 



(^-m + l)!2'"-i(m- 1)! (kr)' 



= (-1) 2 ^ &m = for odd m 

= U, = ^ for even m 

Individual terms of the sum in Equation (jH} diverge as fcr — > 0, but the entire sum 
converges as 



For the first term in Equation (@)), we need the integral 

1 /•°° ,2 k'^ + A . , 

= It' + 

with 

(+) ^ 1 (^ + ^ -!)!(-< + e-^-(P + A) 

^ 4 - m + l)!2'"-i(m - l)!r™ i-oo [(F + yl)2 + B^] k^-^ 

^ 4^J£-m + l)!2'"-i(m-l)!r™i-oo [(P + A)2 + 52] ^"-2 ' 

The contours used for the evaluation of these integrals are shown in Figures la and 
lb. In both cases, the integrand is vanishingly small on the infinite semi- circular 
parts of the contours. On the real k-axis, and on the semicircle around A; = 0, the 
sum of and I^f^ gives us the convergent integrand we need for Equation Q. 
Thus we conclude that 

= 2Tii X [(sum of residues at poles within the contour of Figure la) — 
(sum of residues at poles within the contour of Figure lb)]. (10) 
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The poles within the contour of Figure la are at the zeroes of {k"^ + A)"^ + with 
positive imaginary parts, i.e they are at ±m + iv, where 



u = 



WA^ + B^-A 



V = 



WW+B^+A 



The poles within the contour of Figure lb are at ±u — iv, and at /c = 0. The final 
result is 
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^2-^2 105' 
(A2 + 52)2 

^2-^2 945 
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A(A2 - 3^2) 10395 



■ etc 



(A2 + 52)2 2r5 (^2 + 52)3 

Here Xi represents the contribution to ^ of the residue at the k = pole. 

The k integral for the second term in Equation is performed in the same way. 
Finally, it is necessary to do the x integration in Equation ((Tj). This must be done 
numerically. However, since the range is finite (0 < x < 1) and the integrand is 
smooth, the integrand can be accurately performed with relatively few points. The 
a;-integrations in the results shown below used Simpson's rule, with an x-interval of 
0.001. An example of ug{r), for v = 0.5 and T = 150MeV", is shown in Figure 2. 
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III. Numerical Solutions of Schrodinger Equation 

The binding energy of the quark and the antiquark can be calculated numerically 
by finding bound-state solutions of the Schrodinger equation. In this case, the non- 
relativistic Schrodinger equation governing the relative motion of the two quarks is 
written in the form 



-— + Uo(r) + U2{r)P2{cose)]^lj{r,e,<i)) = E ^{r,e,(j)) 



(12) 



where /i = mQ/2 is the reduced QQ mass. We have included only 'Uo(r) and U2{r) 
in Equation (|T2|l because, in the important r-region, ui{r) with £ > 2 are negligibly 
small (see Figure 2). Because the potential is axially-symmetric, the eigenf unctions 
will be characterized by a definite m— value. However, the spherical symmetry of 
the potential is destroyed by the U2{r)P2{cos6) term, and so the eigenfunctions will 
not be characterized by a unique value of the total angular momentum. Thus a 
solution must be constructed as a linear combination of total-angular-momentum 
eigenstates: 

Mr) ^Hia A\ 



(13) 



If this expansion is substituted into Equation (fT^ the result will be a set of coupled 
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second-order differential equations for the radial functions (f)e{r): 



^-^Uo{r) - k']Mr) - Ti\RMr) M,,, Mr) = 0. (14) 



In Equation (fT^ we have introduced the notation 
E = 

iyL{OA)rY^{e,ct>)Yi{eA) 
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If we use explicit expressions for the vector-coupling coefficients, we obtain 
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In order to describe bound states of the meson we must find normalizeable 
solutions of the coupled equations which are regular at r = 0. This defines 
an eigenvalue condition for k"^. A convenient numerical approach to this problem 
is suggested by the simple special case in which U2{r) = 0. The coupling terms in 
Equation (fT^ vanish, leading to the single equation 



dr'^ 



2/i 



Mo(r) - k^ ] Mr) 



0. 



(15) 



We choose an arbitrary radius, R, and we guess a value for k"^. We then find an 
mtenor solution 0^(r) by numerically integrating Equation (fT3j) from r = to r = R, 
starting at r = with the behavior 

Mr) 



(16) 



Then we find an exterior solution 0^(r) by numerically integrating Equation (fTH|) 
from a very large value of r down to r = i?, starting at the large value of r with the 
behavior 

Mr) ™ h]{zkr). (17) 

The eigenvalue condition on k is that the relative normalizations of the interior and 
exterior functions can be chosen so that there is no discontinuity in value and slope 
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at r = R. This requires that their logarithmic derivatives at r = i? be equal, which 
we can expressed as the condition 



0. 



k"^ in Equation ()15|) is varied until this condition is satisfied. The values of k"^ so 
obtained are independent of the arbitrarily chosen R. This can be seen by using 
(fT^ to show that 
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' dr 
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so that — is independent of r, and so if (fTHj) is true at one value of r, 

it is true at all r. 



To generalize this procedure to the full set of coupled equations (HD), we define 
sets of interior and exterior functions by 



6e^e2 h\{ikr) 



(20) 
(21) 



We now attempt to choose linear combinations of these interior and exterior func- 
tions 



(22) 
(23) 



in order to achieve continuity of value and derivative at the matching radius r = R. 
This requires that the coefficients and satisfy 



(24) 
(25) 



We can express this as a condition on the (5i alone by using Equation ()24|) to eliminate 
from Equation (j^ : 



"^1 = H [ « ^(^) ]h/2 



Y [ h'-\R) a\R) a-\R) h{R) = (3^,. 



(26) 
(27) 
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The necessary and sufficient condition for a non-trivial solution to Equation (|27|) is 

det [ b'-\R) a\R) a-\R) b{R) - 1 ] = 0. (28) 

This is the multi-channel generalization of Equation (jl8|) . The value of fc^ used in 
the numerical integration of the coupled equations (fT^ is varied until Equation (j^Hjl 
is satisfied to an acceptable accuracy. Then Equations and (P7|) determine the 
ae and (3e up to an overall multiplicative factor, which can be obtained from the 
normalization condition 

/•R roo 

E / dr[j:aeMr)a,,r + E / dr [ Y.be,e,{r)(3e, ? = 1. (29) 

As in the one-channel case, the consistency of the procedure guarantees that the 
calculated energy eigenvalues and eigenfunctions are independent of the choice of 
the matching radius. 

The numerical integration of the coupled differential equations was performed 
using the 4th-order Runge-Kutta method with a step-length of 0.01 fm, starting at 
a minimum radius of 0.001 fm [4]. Although the £ sum in Equation ()14p has, in prin- 
ciple, no upper limit, the sum must be truncated in order to make the calculation 
finite. In the results to be shown below, the upper limit on £ was taken to be 8. In 
fact, the strength of the coupling potential U2{r) is small enough so that there was 
very little mixing of £ 7^ components into predominantly £ = eigenfunctions. 

The binding energies obtained in this way are plotted in Figure 3, as functions of 
V and T. Here we have used a heavy-quark mass of irtq = A.SGeV/c^. It is seen that 
our calculations indicate that bound mesons cannot exist at temperatures greater 
than about T = 275 MeV. Furthermore, at any given temperature, the binding en- 
ergy is a decreasing function of v. 
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IV. Dissociation Time of a Heavy Meson Due To Collision With Thermal 
Gluons 

Let us introduce the calculation of the dissociation time by asking: What is the 
rate of photo-ionization if we place a hydrogen atom in a cavity of thermalized 
photons at temperature Tl The Planck distribution [6] states that the flux of 
photons in the interval dk around k is 

^^^^^ 

Let a{k) be the cross-section for photo-emission when the incident photon has mo- 
mentum k. The rate for this process is then 



R 



ck'^a{k) 



:dk. 



/o TT-^le^'^'^/T - 1] 
According to Gasiorowicz, Quantum Physics, Chapter 24 [7] 

Vpee^ 



(31) 



a{k)^ 



2Trh^mc'^k 



(32) 



'4'i{f) is the normalized wave function of the initial bound electron state, Pe is the 
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asymptotic momentum of the outgoing electron [8] , and m is its mass. 



_\ r— >oo 



-ikeV 



Ake-r 



(33) 



~* Pe 

ke — ^ — ke[sm6\cos(j)^ X + sin^'y) + cosO^z]. 
The dfl integration is over 9', 0', the asymptotic direction of the outgoing electron. 

Let us now consider the transition from QED to QCD. In QED the electron- 
photon vertex [9] is associated with matrix clement of 



/47re2 
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(34) 



In QCD the quark-gluon vertex [9] is associated with matrix element of 



(35) 



The A" is a generator of color SU(3). (a =1,2,.. .,8). We are going from a meson 
(color singlet) to a QQ octet. 



< [octet] a I A" I [singlet] >= 



Thus 



which gives us 



3' 



2 

he 3 



The photon flux density is multiplied by 8 to get the total gluon density (8 types of 
gluons, 8 generators of SU(3)). Finally, 



R- r^^^^^dk 

Jo n^[e^kc/T_if''^ 

27rh^ljLek 



and 



(36) 



(37) 



where ji is the reduced QQ mass. 



In his evaluation of a{k), B. MuUer [10] replaced e^^^ by 1 ( long-wave-length ap- 
proximation) and ■0/(r) by e^^'^''^ (no distortion of the outgoing wave) and neglected 
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the connection between kg and k. We will make none of these approximations. 
Now let us expand the final outgoing quark wavefunction 

V V f m ' 



t-pe''''i)Ar) = -e-Ve'^'i)i{r) (39) 

The last three equalities are satisfied since e ■ z = (gluon has transverse polar- 
ization) and since ipi{r) ^ ipi{r) (nearly spherically symmetric ground state). The 
Wi{r) are calculated using the methods of Section III. 

Let 



I{kQ) = J d'rrji^e ■ pe'''^,{r) , (40) 
which is the same as 

I{kQ) = -^Y.^-'e^''yL{d\ct>') f d'rYi\e,<P)^i-re^''^i{r), (41) 
I \/V T:t J r 



where 



e ■ f = e^sinOcoscj) + ej,sin6'sin0 ~ y ~ "^^yWii^, (p) + {ex + iey)Y\( 



Now we use the relation 



-je{kr) 



Yl,{9, <P)e^'^ = ^- ^ r'-'^m + 1)(£ + IY-^YUO, 4>) . (42) 
Substituting Equation (ji^ in Equation (jH}, we get 

47r3/2fc Lmax 

I{kQ) = -=- Y: e^'V£(2£ + l)(£+l)x (43) 
VVk 

[(e, - tey)Y,\e', 0') - (e, + tey)Y\{e', 0')] / drUkr)We{r)^',{r) . 

When we integrate over 9' and cj)' (the dVt integration of Equation (37)), the or- 
thonormality of the spherical harmonics converts the coherent sum over £ in /(^q) 
into an incoherent sum over £ in cr(A;) 
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/' CXJ 

^(2^ + + 1) / drjt{kr)We,{r)ij'-Xr 



(44) 



Equation (j44|l is inserted into Equation (j36|l . and the integration over A; is carried 
out numerically. Figure 4 shows a plot of calculated lifetime r = (1/i?) as a function 
of T and v. 
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The most significant feature of Figure 4 is the prediction of a minimum lifetime. 
When T decreases below about T ~ 170 MeV, the disintegration decreases because 
of falling gluon flux. As T increases above T ~ 170 MeV, the disintegration rate 
decreases because of falling cr(fc). The maximum disintegration rate at T ~ 170 
MeV corresponds to a lifetime of ~ 27.5 fm/c. 

The decrease of cr(A;), in the important k region, with increasing T, is due to 
requirements of conservation of energy and momentum in a process in which a gluon 
is absorbed and a quark is emitted. For larger fc, it is increasingly difficult to find, in 
the meson inital state, quark momenta that are large enough to satisfy both energy 
and momentum conservation. An analagous phonomenon occurs in the photoelec- 
tric effect, and results in photoelectric cross-sections that decrease with increasing 
photon momentum, once a threshold has been passed. In our situation, this effect is 
exacerbated by the decrease of meson binding energy with increasing temperature 
(Figure 3). This results in a meson wave function with greater spatial extent, and 
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thus a lower probability for the occurrence of the necessary high momentum com- 
ponents in the quark-antiquark wave function. 

No lifetime minimum with respect to T was predicted by MuUer [10] since he did 
not take into account both momentum and energy conservation. The expected life- 
time of the quark gluon plasma is 10-20 fm/c, which is appreciably shorter than our 
minimum predicted meson-disintegration lifetime of 27.5 fm/c. Thus our conclusion 
is that the plasma will cease to exist before a significant fraction of the mesons will 
have been disintegrated by gluon impact. Therefore, this process will not have an 
appreciable effect on the number of mesons observed. 

V. Relation of Results to Experiment 

Recent experiments [11] at the CERN SPS have shown a large suppression of 
J/ip production in central Pb+Pb collisions. Following the original idea of Matsui 
and Satz [12] that J/ip would be dissociated in a quark-gluon plasma due to color 
screening, the observed J/%Ij suppression has been suggested as an evidence for the 
formation of the quark-gluon plasma in these collisions [13-15]. 

Since the T meson states in a quark-gluon plasma are also sensitive to color 
screening [12, 16, 17], the study of the T meson suppression in high energy heavy 
ion collisions can be used as a signature for the quark-gluon plasma as well. Because 
the binding energy of T is larger than that of J/iJ, the critical energy density at 
which an upsilon meson is dissociated in the quark-gluon plasma is also higher [18]. 
One therefore expects to see the effects of the quark-gluon plasma on the production 
of the T meson only in ultra-relativistic heavy ion collisions such as at the RHIC 
and the LHC. As in the case of J/ip, one needs to understand the effects of the T 
meson absorption in hadronic matter in order to use its suppression as a signal for 
the quark-gluon plasma in heavy ion collisions. 

It is shown in Figure 3 that a bound T meson cannot exist above a temperature 
of about 275 McV, no matter what velocity it has relative to the plasma. Thus 
observation of T mesons implies that the > 275 MeV plasma never existed, or that 
the mesons were created in a plasma-free region. Bound T mesons can exist at 
lower temperatures, even at speeds relative to the plasma that arc an appreciable 
fraction of c. Thus the complete absence of T mesons in a situation in which they 
should have been created is an indication of the existence of a T > 275 MeV plasma. 

In these considerations, it is not necessary to be concerned about dissociation of 
an T meson due to absorption of a gluon, in a process analagous to the photoelectric 
effect or the photo-disintegration of the deutron. Our calculations show that the 
lifetime of this process is longer than the expected lifetime of the plasma, so once a 
bound T is formed, it has a high probability of surviving the collision. 
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Figure Captions 

Figure 1. Contours for the evaluation of 
Figure 2. ue{r) for T = 150MeV and v = 0.5 

Figure 3. Calculated meson binding energies as a function of T, for various values 
of V. 

Figure 4. Calculated meson dissociation lifetimes as a function of T, for various 
values of v. 
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